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Segment companding laws are being considered for use in PCM to 
cater to the possibility of digital processing of uniformly quantized signals. 
A unified formulation of the segment companding laivs is presented, 
which permits the detailed structure of quantization to be explicitly charac- 
terized. Most important, it is shown that this formulation yields algorithms 
for the systematic synthesis of coders, decoders, and digital compandors. 
Examples of circuits synthesized from the algorithms are shown for two 
segment law companding families — the fi-laiv and the A-law. 

I. INTRODUCTION 

Instantaneous companding of PCM is used to maintain a reasonably 
constant signal to quantizing noise ratio for a wide dynamic range of the 
input voice signal. A logarithmic "/i-law" is a typical "smooth" input/ 
output transfer characteristic that has been approximated in the Dl 
channel bank realized by a diode compandor. 1,2 In contrast to the 
smooth law, segment type companding laws, which are piecewise linear 
approximations to the continuous laws, have also been considered. In 
early experimental PCM systems a segment law was used primarily 
to avoid implementation difficulties associated with tracking of diode 
compandors. 3,4 However, since the mid-sixties segment laws have been 
reconsidered from the view point of digital linearization to cater to the 
possibility of digital processing in an integrated digital network. Two of 
the segment law families have been stressed — the 13-segment A-law, 8,6 
and the 15-segment ju-law used in the D2 channel bank. 7 

The major advantage of the segment laws resides in the digital 
linearization feature, through which digital processing such as digital 
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filtering, conferencing, net loss adjustment, echo suppression, equaliza- 
tion, companding law conversion, and so on, can be performed effi- 
ciently. 8 Also, it has been argued that the segment laws, being well de- 
fined laws, should readily permit an international interconnection on a 
digital basis. However, the defining features of the segment laws have 
not been clearly drawn. 

Many attempts have been made to characterize a particular law 
by table-look-up, by a set of equations, or by forming algorithms, 5,8 
but important parameters defining the particular law have not been 
delineated. What is most important, the lack of an explicit expression 
for the family of segment laws has produced difficulty in the under- 
standing and design of codecs and digital compandors. 

In this paper a unified formulation is presented, which leads to the 
systematic synthesis of codecs and digital compandors for the binary 
code.* In the second section, starting from conventional concepts of the 
/x-law and the A-law, we derive expressions for the decoder output level. 
Output levels are given in terms of a standard form with segment edge 
parameter a and centering parameter c. Coder decision levels may also 
be expressed by the same formulation with the addition of the coding 
parameter 6, thus completely specifying the companding law. This 
approach avoids table-look-up. Relations between the step size, effec- 
tive coder input, transfer gain, and other parameters of the law are 
derived. 

In the third section, we show that this formulation is directly 
related to the codec circuit implementation, and these circuits can be 
synthesized in a systematic and straightforward manner. Unipolar de- 
coders are first synthesized from the formulation, and the synthesis is 
extended to bipolar decoders including the sign bit. We also show 
that the coding parameter b is related to the design of logic in a se- 
quential comparison coder. 

Digital compandors can also be synthesized by using the same formu- 
lation. In the fourth section we display digital expansion and com- 
pression algorithms, and compandor circuits synthesized therefrom. 
At Bell Telephone Laboratories, pioneering studies on digital com- 
pandors have been made by T. P. Stanley, M. R. Aaron and D. G. 
Messerschmitt in unpublished memoranda. W. L. Montgomery has also 
done work in this area. 8 The synthesis technique developed here serves 
to unify these separate works, and makes it possible to design the 
circuits in a systematic way. 

* Extension to arbitrary bases follows in a straightforward manner and will 
not be considered. 
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II. FORMULATION OF SEGMENT LAWS 



2.1 Definition oj Variables 

Suppose that a digital compressed code X is composed of m binary 
digits called "characteristic bits" representing the segment number L, 
and n binary digits called "mantissa bits" representing the quantizing 
step V in a segment, as shown in Fig. 1. The total number M of seg- 
ments in one polarity is equal to 2 m , and the total number N of quantiz- 
ing steps within a segment is equal to 2". Then the digital representation 
of the compressed signal X is expressed in terms of L and V, as 



X{L, V) = V + N-L 



(1) 



or 



V = X mod N 




1 (N-1) 1 (N-1) i....(N-1) 

STEP NUMBER V IN SEGMENT 



I L 



1 N N + 1 2N 2N+1 3N-. 

COMPRESSED SIGNAL X 



Fig. 1 — Segment companding laws — standard form. 
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where X, L and V are denned in the regions as 

X = digital compressed signal £{0,1, • • • , MN — 1}, 

L = segment number £{0, 1, • • • , M — 1}, and ? (2) 

V = step number within a segment £ {0, 1, • • • , N — 1}.J 

The definition of X, L and V beginning from "0" permits a consistent 
algebraic treatment and yields a straightforward circuit design. 

If we denote the decoder output level by y , then it is a function of 
the digital compressed signal and y = Vo(X) = y (L, V). The quantiz- 
ing step size is then given by the difference of two adjacent output levels, 
and is 

&o(L) m y (L, V + 1) - y (L, V), V * N - ll (g) 

K(L) = y (L +1,0)- y (L, N - 1), V = N - lj 

where A' Q (L) is the step size at the segment edge. Within a segment the 
step size A (L) is constant and independent of V. 

2.2 Definition of the Segment ix-law 

We will define the segment ju-law and derive an expression for an 
output level as a function of L and V. The "so-called" segment /z-law 
should satisfy the following: 

(i) Except for the segment edges, the ratio of the step sizes of the 
adjacent segments is equal to 2; that is, 

A ° ( a^ 1] = 2 < £€{0,1, •■■ ,M-2], (4) 

(it) y (0, 0) = c (centering), (5) 

(Hi) A (0) = 1 (normalization) . (6) 

Solving the difference equation (4) subject to equations (3), (5) and (6), 
we have 

Vo (L, V) = 2 L (V + N + a)-N-a + c 

= UL, V + a) - a + c (7) 

for an arbitrary real a, where / M (L, V) is the standard form of the seg- 
ment /i-law, and is 

UL, V) m 2 L (V + N) - N. (8) 

The standard form is interpreted in the following way. If we take a 
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continuous variable v such that 

v £ [0, N] 

then j{L, v) represents a continuous segment /x-Iaw as shown in Fig. 1. 
The continuous segment p-law / M (L, v) is a set of segment lines whose 
segment edges / M (L, N) = /„(L + 1, 0) satisfy the so called /i-law 1 for 
x = LN + v, namely: 

x_ = log (1 + »y/E 2 ) 
Bi log (1 + M ) w 

for normalization constants #, = MN, E 2 = N(2 M — 1), and m = 
2 — 1. In other words, the continuous standard form gives an exact 
"chord type" approximation to the smooth yi-law. On this continuous 
segment law the discrete point V should be defined as the truncation of v. 
The output level (7) is then the modified version of the standard 
form with parameters a and c, and it constitutes a class of "/i-like" 
laws. The centering parameter c determines the origin, and the edge 
effect parameter a describes the discontinuity of step size at the seg- 
ment edges. The continuous segment law varies with parameter a as 
shown in Fig. 2, and produces discontinuities at the segment edges. 

2.3 Extension to the Segment A-law 

To obtain the segment A-law, 5 we replace the condition (i) by 



/ 



/ 
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y// L = 2 






// * 




,/ 
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y v 
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L = 1 



L = o 



COMPRESSED SIGNAL, X 



Fig. 2 — Effect of segment edge parameter a: Curve 1 — a > 1; Curve 2 — a = 
and Curve 3 — a < 0. 
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(i)' Except for the segment edges the ratio of the step size of the 
adjacent segments is 

Ao(L+l) _ [l, L = 0, (10) 

Ao(L) [2, LG {1,2, ...\m -2}. 

This implies that the 0th and 1st segments constitute a long colinear 
segment and the step size ratio of the adjacent segments elsewhere is 
equal to 2. Then solving equation (10) subject to equations (3), (5) 
and (6) we have 

y (L, V) =f A (L,V + a) - a + c (11) 

where f A (L, V) is the standard form of the segment A-law and is 

U(L,V) = \ V > L = °' (12) 

Note that equation (11) is identical with equation (7) for the /t-law, 
and the output level for the A-law is characterized by the standard 
form f A (L, V) in equation (12). 

As in the ^-law case, where a discrete V is replaced by a continuous 
variable v £ [0, N] in the standard form f(L, V), a "chord type" approxi- 
mation to the so called A-law is obtained. The A-law has been de- 
fined as 5,9 



x 

A' 



, f frff> ■ ve[o,E 4 /A} 

1 + log* A 

1 + lo S . (Ay/E,) ye[Ei/A:Et] 

1 -f log A - A 



(13) 



for K = 2, A = 128,* E 3 = MN, and E A = N -2 M . 

Expression (11) applies for any type of segment companding law if 
the standard form j A (L, v) is appropriately defined to yield a continuous 
piece-wise linear approximation to a continuous curve. 

2.4 Decision Levels 

In the foregoing we have concentrated on the decoder output level 
because the level of primary importance in companding is the output 
level. The coder decision level should then be determined to yield an 
appropriate range of the input signal corresponding to that output level. 

* The parameter values A = 87.6, K = e (base of natural logarithms) are often 
used. 6 However, to be consistent with a "chord type" approximation, K = 2 and 
A = 128 are used here. 
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First, the decoder output level should be a good representative of the 
signal range between two adjacent decision levels. It has been shown 
elsewhere that for a given set of output levels the locally optimum 
decision levels, which minimize the mean squared error, should lie 
halfway between two adjacent output levels irrespective of the input 
signal distribution. 10 This can be achieved, except at the segment 
edges, by shifting the output level a half step up or down in the com- 
pressed signal domain; that is, the decision level is 

y d (L, V) = y (L, V + b) (14) 

where b is the coding parameter and takes on the values +0.5 or —0.5. 
When b — +0.5, y d (L, V) gives the upper bound of the input range 
corresponding to X(L, V); when 6 — —0.5, it corresponds to the 
lower bound. The coder decision level y d implies that an analog input y 
is encoded into X(L, V) if 

y d (X - 1 | 5 - +0.5) =g y < y d (X \ b = +0.5) (15) 

for the upper bound approach (6 = +0.5), and if 

y d {X | b - -0.5) ^ y < y d (X + 1 | b = -0.5) (16) 

for the lower bound approach (b = —0.5). In general, these two ap- 
proaches are identical at all points except at segment edges. Now, for 
any points including segment edges, it can be shown (Appendix A) 
that the two approaches coincide if and only if the edge effect parameter 
a is equal to 0.5. As we will see in Section 3.4, the coding parameter 
b is closely related to the binary coding logic in a sequential comparison 
coder. 

2.5 Codec Transfer Characteristics 

The decoder output level and the coder decision level have been 
denned. Now, the characteristic from the coder input to the decoder 
output as shown in Fig. 3 can be expressed explicitly by two equations 
with the auxiliary variables L, V and the parameters a, b and c; that is, 

y (L, V) = 1(L, V + a) - a + c 1 (1?) 

y d (L, V) = f(L, V + a + b) - a + c\ 

where /(L, V) is the standard form defined by equation (S) for the 
M-law and equation (12) for the A-law. These equations characterize 
the precise quantization structures. 

The centering parameter c is or 0.5, and when it is the transfer 
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characteristic at the origin is "mid-tread" as shown in Fig. 3a. When 
c = 0.5, it is "mid-riser" as shown in Fig. 3b. 

To show the effect of the parameter a, a few examples of the transfer 
curves are shown in Fig. 4 for N = 4 and mid-tread operation (c = 0). 
As mentioned before, when a = 0.5, the transfer curves for b = +0.5 
and —0.5 are identical except for a shift in the numbering of X(L, V). 

2.6 Step Sizes 

The step sizes of the output levels and the decision levels are given, 
except for the segment edges, by 

Ao(L) = y (L, V + 1) - y (L, V)] y ^ N _ ± (18) 

A d (L) = y d (L, 7+1)- y d (L, V)\ 
and the corresponding step sizes at the segment edges are given by 
K(L) = y (L +1,0)- y (L, N - 1),| (1Q) 

A' d (L) = y d (L +1,0)- y d (L, N - 1).J 
It is easily shown that for non-edge points 

AoOO = A d (L) = 2 L (20) 

for the ju-law. However, for the edge points the step size is given by 

A> Q (L) = /(L + 1, a) - f(L, N - 1 + a) = A (L)-(1 + a),' 

A' d (L) = f(L + 1, a + 6) - f(L, N - 1 + a + b) 

- A d (L)-(l + a+ 6).. 
Step sizes at the segment edges are depicted in Fig. 5 with respect to a. 



(21) 



D 
n 
h- 

8 1 



T7 



(a) 



Z 



z 



0*1 2 3 

INPUT y 



3 


(b) 


A 


U) 




/ 


1- * 

D 

a. 

H 




f 




O ' 


/ 






7^ , , 


j 





I 

NPUT 


l 3 

y 



Fig. 3— Effect of centering parameter c, (a) mid-tread (c = 0), and (b) mid-riser 
(c = 0.5). 
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Fig. 4 — Transfer staircases with N = 4 and mid-tread for (a) RLA (a = 0, b = 0.5), 
(b) DLA (a = 0.5, & = 0.5), and (c) RLA (a = 1, b = -0.5). 



A (L) = A d (L) - 



For the A-law, equation (20) is modified to 

1, L = 0, 

and relation (21) applies for the A-law on the understanding that the 
"segment edges" exclude the colinear segment edge between 0th and 1st. 

2.7 Effective Coder Input 

Now, we will define an effective coder input y(L, V) such that y(L, V) 
is the statistical average of the input signal between two adjacent 
decision levels. If we assume that the input is uniformly distributed 
between two adjacent decision levels, then y(L, V) is equal to their 
average as shown below. For b = -f-0.5 (upper bound case) 

V)], V 5*0, 

ih[y d (L - 1, N ■ I 1- ya(L, 0)], 7 = (edge). 



1564 THE BELL SYSTEM TECHNICAL JOURNAL, SEPTEMBER 1970 



2.5 



UJ 1.0 



2 
^^ 3 



0.5 1.0 

SEGMENT EDGE PARAMETER, <2 



Fig. 5 — Step size ratio at the segment edges relative to the lower segment step 
size: Curve 1— output level; Curve 2— decision level {b = 0.5); and Curve 3— 
decision level (b = —0.5). 



For b = — 0.5 (lower bound case) 

\[y d (L,N-l) + y d (L+l,G)}, V = N - 1 (edge). 

Except for the edge points, we have 

y(L, V) = y (L, V) (22) 

and this is a direct consequence of our definition of y d in equation (14). 
At the segment edges (L„ , V „) , it can be shown that 



y(L. , V.) = y a (L, , V.) - A„(L)-6-(a - 0.5) 

where edge point means 

[(X + 1, 0) when 6 = +0.5, 
\(L, N - 1) when b = -0.5. 



(L t , V.) = 



(23) 



(24) 



Therefore, except for the segment edges the transfer gain from the 
effective coder input to the decoder output is unity. At these segment 
edges, if we define the tracking error e, which is the deviation from 
unity of the transfer gain from the effective coder input y to the decoder 
output 7/ n , then e is 
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e = 1 — y/yo 

(An(L) 



= 1 2/o 



b-(a — 0.5), at segment edges 



(25) 



0, 



elsewhere. 



The ratio e/(Ao/z/o) is the gain error relative to the step size, and is 
depicted in Fig. 6, where it is soon that the tracking error is zero over 
the entire input range if and only if a = 0.5. 

2.8 Relation Between Parameters 

In Section 2.4 it was noted that for a given output level the locally 
optimal (minimum mean squared error) decision level is half way be- 
tween adjacent output levels. 10 It can be shown that this optimality 
condition is satisfied for any X including segment edges, if and only 
if a + b = 0.5, which means a = for b = 0.5, or a = 1 for b = —0.5. 
Since optimal placement is made from the given output level, this is 
called "Reconstruction (or Output) Level Assignment" (RLA). 8 

On the other hand, for given decision levels, assuming the uniform 
distribution of the input in each quantization step, the locally optimal 
(minimum mean squared error) placement of the output levels lies 
halfway between adjacent decision levels. 10 It can also be proven that 
this condition is satisfied for any point X including segment edges, if 
and only if a = 0.5, irrespective of b = ±0.5. Since the optimal ap- 
proach is based on the given decision levels, this is called "Decision 
Level Assignment" (DLA). 8 Proof of the above statements is given in 
Appendix B. 






0.5 1.0 

SEGMENT EDGE PARAMETER, Q. 



Fig. 6 — Tracking error at the segment edges. 
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As shown previously, when a = 0.5 (DLA) the effective tracking 
error is zero for all inputs, and the upper and the lower bound ap- 
proaches give identical coding. Furthermore, as will be seen later, the 
digital compression algorithm and circuits are simpler in DLA than 
in RLA. The D2 channel bank has been designed according to DLA. 
Hereafter, our attention is mainly focused on DLA, although the 
analysis is easily generalized to include other cases. 

When M = 8, the segment /x-law is generally referred to as the 
"15-segment n-\&w" including negative segments. In contrast, the 
segment A-law with M = 8 is called the "13-segment A-law." In 
Table I the output levels and the step sizes are listed for the 8-bit 
(including sign bit) 15-segment ju-law and 13-segment A-law. The 
maximum and minimum output levels are shown in Table II for ready 

reference. 

Quantizing noise, as calculated for the Laplacian (negative expo- 
nential) distributed speech input y in ,yields the output signal to quantiz- 
ing noise ratio shown in Fig. 7. Signal and noise are denned in a band 
of width equal to half the sampling rate. The average transfer gain (G) 
is chosen so as to minimize the quantizing noise, and is given by 
E[y in -yo]/E[y 2 in ]. The idle channel noise, calculated for a zero-mean 
gaussian input noise, is simply ICN = E[if ] and is shown in Fig. 7. 
The centering parameter c is important in the small signal range as is 



Table I — Output Level Step Sizes for m = 3 and n = 4 
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Table II — Maximum and Minimum of Output Levels for m = 3 

AND 71 = 4 
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Fig. 7 — Average gain, signal-to-noise ratio, and idle channel noise for DLA with 
m = 3 and n = 4: Curve 1 — /i-law, mid-tread (c = 0); Curve 2 — /i-law, mid-riser 
(c = 0.5); Curve 3 — A-law, mid-tread (c == 0); and Curve 4 — mid-riser (c = 0.5). 
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well known, and it affects center clipping and idle channel noise. The 
effect of a on the signal quality is very small. The signal-to-noise ratio 
relative to the S/N for DLA case (curve 1 in Fig. 7) and the average 
gain are shown in Figs. 8a and b as curves 1 and 2 for the two RLA 
cases. The small difference (±0.14 dB) in S/N in the lower signal range 
is due to the difference in overload levels. However, if there is a mis- 
tracking such that the signal is encoded with a, = 0.5 and decoded 
with a 2 = 1, the mismatch results in distortion and net loss variation. 
The penalty in the signal-to-noise ratio and the change in average gain 
are shown as curve 3 in Figs. Sa and b. Although the penalties are 
very small, the net loss variation barely meets toll quality requirement, 
and accumulated impairment in multi-link digital conversions would 
be intolerable. Therefore, the mixed use of DLA and RLA is to be 
avoided. 

2.9 General Formulation 

To summarize this section we give a general expression which is 
conveniently used hereafter in the synthesis of codecs and digital 
compandors. From equations (7), (8), (12) and (14) the general ex- 
pression can be written as 



y{L, V) = A(L)-(7 + P) - Q 
where for the segment M-laws 
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Fig. 8 — Effects of parameters on S/N and average gain for /z-law with m = 3 and 
n = 4: Curve 1— RLA (a = 0, b = 0.5); Curve 2— RLA (a = 1, b = -0.5); and 
Curve 3 — mistracking case coded by DLA and decoded by RLA (a = 1). 
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(27) 



Equation (26) can represent either the output level by equating 6 to 
or the decision level with b = ±0.5. For the segment A-laws, expression 
(27) is replaced by 



(28) 



where 



A(L) = L 


- V, 


P = N- 


■q + a + b, 


Q = a- 


~ c, 




L = 0, 
L ^ 0. 



III. SYNTHESIS OF CODECS 

3.1 Synthesis of Unipolar Decoder 

In this section, we show that the representation of the output level 
and the decision level by equation (26) leads directly to the implementa- 
tion of the decoder. In the following, attention will be focused mainly 
on the segment p-law. Extension to the A-law is readily accomplished 
by similar procedures. 

Equation (26) given by 

y(L,V) = A(L)-(7 + P) -Q 

is clearly interpreted as adding bias P to V, amplifying it by A(L), 
and subtracting bias Q from this result, thus obtaining the decoder 
output y(L, V). 

Let us represent the nonlinear code X by an (m + n) bit binary 
sequence (except for the sign bit e ), e, £ {0, 1}, 

X = X{ei , e 2 , • • • , e m , e m+1 , • • • , e m+n \ 

= £ 2 m+ "- < .e i . (29) 
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Then, 



L = L{e x ,e 2 , ••• , e m ) = £ 2 m 'e, , 

(30) 

V = 7{e m+I , ••• ,e m+n ] = Er^.- • 

i-l J 

From equation (27), A(L) can be written as 

A(L) = 2 L = ft {2 (2m " , ' > r < ( 3l ) 

t-i 

and is realized by a cascade of switched amplifiers (or attenuators) 
whose gain is switched between 2 <2m " i) and 1 according to the digit 
e< = 1 or 0, respectively. Likewise, V in equation (30) is obtained 
directly from a uniform decoder of any of the known forms 11 operated 

by {e m+ i , e m+2 , • • • , e m+n }. 

An example of the /x-law decoder is shown in Fig. 9a. Biases P and 
Q are given by (N + a + b) and (N + a - c), respectively. The decoder 
output is either the decoder output level y (L, V) by replacing 6 = 
in P, or the coder decision level y d (L, V) by choosing b = ±0.5. 

The decoder circuit for the segment A-law is synthesized in the same 
way using equation (28), and is shown in Fig. 9b. 

3.2 Consideration of the Sign Bit 

In the preceding, the sign bit e has been excluded from considera- 
tion. In practical cases it is desirable to use a bipolar decoder to cover 
the whole positive and negative signal ranges. We will give alternative 
representations of the decoder output in terms of the signed input 
codes. First, assume that the output level is symmetric with respect 
to the 0-level line.* Suppose that the sign bit e is equal to for the 
positive signal region, and 1 for negative region. If we define 

«o= 1 -2e G {+1,-1}, O 2 ) 

then from equation (26) the bipolar decoder output is given by 

y(L, V) = a {A(L)-(V + P)-Q} (33) 

= A(L){a V + a P] -a Q. (34) 

If we look at the binary representation of a V = {e , e m+i , e m+2 , 
••• , e m+n ) as in equation (30), the so-called "sign and magnitude" 

* This implies that in the mid-tread case the "0" output level is shared between 
the positive "0" and the negative "0". 
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Fig. 9 — Synthesis of segment law decoders (m = 3, n = 4): (a) 15-segment /x-law 
and (b) 13-segment A-law. 

or "folded binary" representation results. Another representation often 
used in computer arithmetic is the "l'fl complement" code. 12 If we 
denote the l's complement code by 

V = V'{e' mli ,e' m+2 , ••• ,e' m+n ] = Z 2»-*< +i , 



then V is related to positive V by 

a V = V - e n (N - 1) 
and also 

Oi = 6 m+ , e mod 2. 



(35) 
(36) 
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Therefore, equation (34) can be expressed in terms of V and 

y(L, V) = A(L){V - e (N - 1) + a n P} - a Q. (37) 

Note that the modification in the representation applies only to the 
mantissa bits V, and the characteristics bits L should remain unchanged. 

3.3 Synthesis of Bipolar Decoder 

We can synthesize the bipolar decoder circuits according to the 
two alternative representations (33) and (37).* 

(i) Sign-and-magnitude (folded binary)— The "sign and magnitude" 
decoder according to equation (35) is obtained by inverting the decoder 
output of the unipolar decoder, or by inverting all the supply voltages 
to the unipolar decoder as shown in Fig. 10a. The former requires a 
precise analog inverter. In the latter circuit, the source switch must 
be in serial with the decoder switch, and stringent requirements arise 
due to the wide dynamic range of these analog switches. 

(m) l's complement — Suppose that we use a bipolar source (+#, — E) 
instead of unipolar source. (+E, 0) By converting the unipolar code 
to a bipolar code, we have 

«Cn = 2< +i - 1G {-1, +1}, (38) 

and also using equation (32), equation (37) becomes 



- a Q (39) 



y(L, V) = A(L)[|: 2 n ~ i - x -a' m+i + a (l> + ^y^). 

where A(L), P and Q are given by equation (27). This leads to the 
decoder circuit shown in Fig. 10b. The weighting currents of the decoder 
are given by the coefficients of the a's in equation (39). The weighted 
current [P + (N — l)/2] that is switched by a is about six times 
greater than the most significant digit weight (= 2 n-2 ) of the linear 
decoder. Therefore, the problem of the serial switch in Fig. 10a is 
replaced by a threefold increase in the accuracy requirements of the 
weighting network here. This appears to be a worthwhile trade in 
the practical design of codecs. 

3.4 Interpretation of Coding Parameter b 

As discussed in Section 2.4, the coder decision level y d is related 
to the decoder output level y by 

V*(L, V) = y (L, V + b) 



* We have derived another representation for the "2's complement code" but 
it is omitted because of its increased complexity over the other two representations. 



SIGN BIT e r 



SEGMENT COMPANDING 

LINEAR DECODER 

± 1 /o ±2 /o ±4 /o i8 /o 



1573 




I 1 i 1 1 



FOLDED BINARY CODE 



(a) 



LINEAR DECODER 
+ °- 5 /-0.5 +1 /"1 +2 /"2 +4 /" 




SOURCE 



I'S COMPLEMENT CODE 



CbJ 
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for b = ±0.5. When b = +0.5 the decision level is an upper-bound 
for the corresponding output level, and when 6 = —0.5 it is a lower 
bound, in the sense given by inequalities (15) and (16). 

In the following we show how the lower or upper bound approach 
is related through the parameter b to the logic used in implementing 
a sequential comparison coder. 

Referring to Fig. 11, for the lower bound approach (6 = —0.5) 
(from the "all reset" state "000"), the first digit trial "1" is made; 
and the decision level corresponding to "100" is compared with the 
input. As a result, if the error e defined as 



i = input — local decoded signal 
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Fig. 11 — Upper bound and lower bound decision levels. 

is negative the trial "1" is "reset" to "0", and remains unchanged 
if e is positive. This procedure is repeated digit by digit from the MSD 
to LSD, and the coding is performed sequentially so as to minimize e. 
On the other hand, for the upper bound approach (6 = +0.5), 
the "all reset" is made to "1", the "trial" is made to "0", and the 
"reset" is made to "0" if e > 0. Therefore, if we take a package of 
"all reset-trial-conditional reset", it is "0-1-0" for the upper bound 
approach, and "1-0-1" for the lower bound approach. From this 
package the circuit design of sequential coding logic and registers 
follows immediately. 

3.5 Synthesis oj Sequential Comparison Coder 

From the result of Section 3.4, and using the unipolar decoder 
described in Section 3.1, the synthesis of a unipolar coder is straight- 
forward. Here, we treat only the bipolar coder. It has been shown 
in Section 3.3 that the bipolar decoder in l's complement form has 
an advantage in circuit design. If the upper bound approach is taken 
(b = +0.5), the "all reset-trial-conditional reset" package is "0-1-0" 
for a positive signal. For the negative half if b = +0.5 the package 
becomes "1-0-1" because of the l's complement (inverted) representa- 
tion. However, the logic package can be uniformly "0-1-0" for the 
entire range, if we take b = +0.5 for positive inputs and 6 = —0.5 
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for the negative region. Figure 12 shows the schematic of a bipolar 
sequential comparison coder based on the l's complement code. The 
local decoder, which must have the nonsymmetric feature in b, is 
obtained from the l's complement decoder of Fig. 10b with 6 = 0. 
In addition, a constant bias b = +0.5 is added, externally. Then, 
the register and logic are operated on the basis of the l's complement 
code, and the coding is performed with the logic package of "0-1-0." 
The segment characteristic bit should always be a "folded binary 
code", which is simply the inverse of the l's complement code. Of 
course, there are many other alternatives for building a coder, and 
the synthesis of any of those immediately follows from our formulation. 

IV. SYNTHESIS OF DIGITAL COMPANDORS 

In Section III we considered the conversion from analogue signal 
to digital compressed code or the converse. Here, we consider the 
conversion between the compressed code and the linear code. The 
conversion from the compressed code to the linear code is called "digital 
expansion" and the converse is called "digital compression." In Sec- 
tions 4.1 and 4.2 the digital expandor synthesis is treated, and after 
that the digital compressor is studied. 

4.1 Digital Expansion Algorithm 

Suppose that a nonlinear code X(L, V) with bit length (m + n) 
is digitally linearized to F. The linearized code Y should correspond 
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Fig. 12 — Bipolar sequential comparison coder with l's complement code logic. 
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to the decoder output level 

Y = y (L, V) = i(L, V + a) - a + c. (40) 

Suppose that the companding law is symmetric for positive and negative 
signals. We generally exclude the sign bit from consideration. In the 
case of RLA and mid-tread (c = 0), 7 is an integer and is expressed 
by a linear code with length (2 m + n) bits. In the cases of DLA (a = 0.5) 
or the mid-riser (c = 0.5), Y is expressed by (2 m + n + 1) bits including 
one fractional digit. For example the D2 code (a = 0.5, c = 0, in = 3, 
and n = 4) is basically* expressed in terms of 13-bit binary code without 
the sign bit. 

The expanded signal Y should correspond to the output level. There- 
fore, from equations (26) and (27) with b = 0, Y is given by 

F = A(L)-(V + P) - Q (41) 

where A(L) = 2 L is a shifting operator, by which the binary sequence 
(7 + P) is shifted to the left by L bits. Therefore, the algorithm for 
code conversion from X(L, 7) to Y is 

Algorithm 1 (expansion): For given X(L, 7), Y is obtained such that 

(i) Add P to 7, 

(m) Shift (7 + P) by L bits to the left, and 
(Hi) Subtract Q, 

where P = N + a and Q = A^ + a — c. All variables and constants 
are assumed to be in binary representation. 

4.2 Synthesis o] Digital Expandor 

The above algorithm leads directly to a parallel expandor circuit 
as shown in Fig. 13. The input L is written into a binary counter BC, 
and the input 7 is written into a shift register SR through a parallel 
adder with addend (N + a). Then, the binary counter BC counts 
down to "000", and stops. This interval is equal to L. The SR is shifted 
by the same clock pulses, and L shifts to the left is performed. Then, 
the entry of SR is read out through a parallel subtractor with sub- 
trahend (N + a — c), thus yielding the parallel expanded output Y. 

In practice, a digital expandor is generally followed by a digital 
processor, in which serial arithmetic logic is most conveniently used 
when the signal sequence begins with the LSD (Least Significant 
Digit). 12,13 In the sequel, we concentrate on the serial implementation. 

* Except for bit robbing for signalling and framing. 7 
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Fig. 13 — Parallel digital expandor. 

If we replace (X2) by an operator z, which at the same time represents 
the delay operator by one clock interval, then Y is given by 

Y(z) = z L -{V(z) + P(z)} -Q(z) (42) 

where V(z), P{z) and Q(z) are binary sequences beginning from the 
LSD. Therefore, a general form of digital expandor can be described 
schematically as shown in Fig. 14. 

Let us consider more specifically the /x-law expandor design. Using 
the 2-transform notation, equation (42) can be written as 

Y(z) = z L {V(z) + z n + a{z)\ - {z n + a(z) - c(z)} 

= z L \V(z) + z"} + a{z)-(z L - 1) + c(z) - z n (43) 

where (±) implies addition/subtraction with carry, if necessary. As- 
suming DLA [a(z) — i" 1 ], mid-tread (c = 0), m = 3 and n = 4, we 
have 

zY(z) = z L [zV(z) + zP(z)} - zQ(z) 

= z L (z° + e 7 z + ej? + e5 z 3 + e^ + z 5 ) - z° - z 5 . 

This can further be written in the form 

zY(z) = z L (e.z + e z 2 + e 5 z 3 + e iZ A + z 5 ) + (z L - 1) - z\ (44) 

From this a serial expandor circuit follows as shown in Fig. 15. The 
inputs L and V are written into a binary counter BC and a shift register 
SR, respectively. Also "1" is written in the last stage of SR, and the 
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initial entry of SR, which is equal to (N + V), corresponds to the 
first term in equation (44) without delay z L . The binary counter BC 
counts down until it becomes "000" and stops, thus giving a gating 
signal T L with length L. Instead of having a separate variable delay 
as in Fig. 14, SR is designed to wait for L intervals before shifting 
to the right. This gives a delay corresponding to z L , and the output 
of the shift register U, which is equal to the first term in equation (44), 
is depicted as shown in the timing chart of Fig. 15. Next, notice that 
the second term {z L — 1) in equation (44) is a sequence of successive 
l's with length L. It is exactly the same as the gating signal T L . Since 
T L does not overlap with the first term U, it can be simply added 
through an OR circuit to the output of SR. The last term z", which 
is nothing but a timing pulse at the 5th time slot, is subtracted through 
a half subtractor. Thus, the output LSD-first sequence z-Y(z) is ob- 
tained immediately after the write-in of the input. 

4.3 Digital Compression Algorithm 

In this section, we consider the conversion from a linear code Y 
to a nonlinear code X. Suppose that a linear code Y is a positive real 
number representing a signal amplitude normalized by the minimum 
quantizing step of the compressed code. This code Y is not necessarily 
an integer, and in general it has q fractional digits in binary representa- 
tion. This corresponds to the common practice to have more digits 
in a digital processor than in the input/output digits to reduce the 
effect of rounding or truncation in arithmetic units (such as digital 
filters). 

As described in Section 2.4 nonlinear coding from an analog signal 
y to a compressed code X(L, V) is performed by comparing the input 
with the decision levels. Similarly, in equations (15) and (16), on 
replacing the analog input y by a discrete linear code Y, we have the 
following basic digital compression algorithm. 

Algorithm 2 (compression: basic): 

(i) Upper bound approach — Choose X such that 

y d (X - 1 | b = 4-0.5 )g Y < y d (X \ b = +0.5). (45) 

(ii) Lower bound approach — Choose -X" such that 

y d (X | b = -0.5) ^ Y < y d (X + 1 | b = -0.5), (46) 

where y d is exactly the decision level defined by equation (14). 

It has been shown that both approaches yield identical outputs 
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if and only if a = 0.5 (DLA). Also note that there is an equality only 
in the lower boundaries. There is the alternative algorithm in which 
the equality in equations (45) or (46) is given to the upper boundaries. 
As compared with the coding from analog input y described in Sec- 
tion 2.4, Algorithm 2 is equivalent to analog coding with input y = 
Y — 2"*"\ and the alternative upper boundary algorithm is equivalent 
to analog coding with y = Y + 2""" 1 . Both algorithms are asymp- 
totically equivalent to analog coding for large q. We use the former 
since it is more directly related to circuit implementation. 

Though this algorithm is basic, it is not convenient for synthesizing 
the circuit. We will divide the above algorithm into the X-search 
process and the ^-determination processes. Let us introduce a new 
variable W such that 

Wit, Y) = IA(L)}~ 1 '(Y + Q) (47) 

where Q = N + a - c. Then, for the DLA the compression algorithm 
is given by 

Algorithm 3 (compression: DLA): In the case of DLA (a = 0.5), the 
conversion is performed as follows: choose L such that* 

W(L, Y) G [N, 2N) (48) 

determine V such that 

V = 3[W(L, Y)] - N. (49) 

Here, 3[-] is the truncation function and denned as 
3[x] s truncation of x 

= integer I such that I g x < I + 1. (50) 

We also define the rounding function 

(R[x] = rounding of x 

= 5[3 + 0.5]. (51) 

The derivation of this algorithm is given in Appendix C. If we take 
the inverse of 

Y = A(L)(V + P) - Q 
with respect to V, then 

V = [AiLyr-iY -f-Q) - P. 



* Semi-closed interval [N, 2N) means N < W < 2 A'. 



SEGMENT COMPANDING 1581 

The first term is exactly W(L, Y) and P = N + a = N + 0.5 in DLA, 
and 

V = W(L, Y) - N - 0.5. 

By taking the rounding of V and using equation (51), we have 

CR[7] = (R[W(L, Y) - 0.5] - N 

= 3[W(L, Y)] - N. 

This is the same as V in the algorithm. Therefore, the physical implica- 
tion of the algorithm is that V is obtained by simply rounding the 
transjormed variable V. 

For the Reconstruction Level Assignment (a = for b = +0.5, 
or a = 1 for b = —0.5) the algorithm is more complicated, 14 and is 
Algorithm 4 (compression: RLA): For the RLA, conversion is ac- 
complished in the following way: choose L such that 

W(L , Y) G [N, 2N), (52) 

determine V~ such that 

V = 3[PF(Lo , Y)] - N (53) 

and correct by 

X(L, V) = Xo(L Q ,Y )-r-e (54) 

where the corrector is given by 

= W of - a (55) 

for W of , the highest fractional binary digit of W(L , Y). 

The proof is complicated, and so is the circuit, and both are omitted 
because of less interest in RLA. 

4.4 Synthesis of Digital Compressor 

We will synthesize a digital compressor circuit according to the 
above algorithms. For simplicity DLA is assumed. Also assume that 
an input Y is a binary sequence beginning with LSD. In Algorithm 2 
the L-search process by equation (48) is accomplished, first by forming 
a binary sequence (Y 4- Q), and then shifting bit by bit to the right 
until equation (48) is satisfied. The criterion (48) implies that the nth 
digit of W is "1", and the digits of W higher than the nth digit are 
all"0". 

Based on this algorithm an example of the M-law digital compressor 
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for DLA, m = 3 and n = 4 is synthesized as shown in Fig. 16. The 
input Y(z), consisting of 12 integer digits and q fractional digits, is 
fed into a serial adder and Q = (N + a - c) is added. Then, the output 
of the adder is written into the first shift register SI^ and then to the 
second shift register SR 2 . Let us denote the time instant by t Q when 
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the MSD of the input sequence is just stored in the first cell a x . While 
the write-in proceeds, the counter BC, which is being counted down 
by the clock, is reset to the "111" state whenever the entry of the 
first cell a, is 1. Therefore, at the instant t the BC entry is exactly L, 
and it is read out at t . This also means that the highest order "1" 
is located at the cell a 8 -L , and an additional L shifts forces the last 
entry a 8 of SR, to be "1". This is accomplished by shifting SR, and 
SR 2 until the entry of BC becomes "000" at time tt . 

From Algorithm 3 the mantissa output V is given by 3[W] — N. 
Circuit-wise, 3[PP] amounts to discarding the smaller digits below the 
last stage a, 2 of SR 2 , and (— N) means discarding a 6 = 1. Thus, t^ 
is obtained from the cells (a 9 , a in , a„ , a 12 ) of SR 2 , which is being 
stopped after L shifts. The L output has been read out at time t from 
BC, thus completing the conversion. 

For mid-riser formulation, the circuit can be modified by changing 
the addend Q = N -\r a — c. A compressor for RLA is also achieved 
according to Algorithm 4 by adding a correction circuit by 8 after 
the circuit of Fig. 16. The correction circuit may be a parallel adder or 
a 7-bit binary counter, and obviously the RLA compressor is more 
complicated than the DLA one. 

Among many alternative circuits, a notable feature of the circuit 
in Fig. 16 is that it can process the input sequence in "real time." 
While the conversion for a code word is going on, the next code word 
is being shifted into the SR, to prepare for the conversion of the next 
code word. One might argue that if an input to SRi is a sequence 
beginning from MSD the shift is made from the right to the left simply 
to yield equation (48). However, doing so one would need additional 
memory to convert the LSD sequence to the MSD sequence before SR, . 

4.5 Extension to the Segment A-law 

The synthesis procedure can be extended to the segment A-law 
using expression (2S) instead of (27). The expansion algorithm (cor- 
responding to Algorithm 1) for the A-law is obtained by replacing 
P = N -r] + a, Q = a — c and the L-shif ts by (L — tj) shifts, where 
t) is equal to for L = and 1 otherwise. Following the same procedure, 
a digital expandor circuit can be synthesized as shown in Fig. 17, 
for DLA, mid-riser, m = 3 and n = 4. The A-law expandor differs 
from the p-law expandor in Fig. 15 in several respects. The input to 
the first cell is "tj", which is obtained from an OR circuit operating 
on e, , e 2 , and e 3 . The input to the last cell of SR is "1". The gating 
pulse T L - n with length (L — ??) is obtained by OR operation on e x 
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Fig. 17 — Thirteen-segment A-law digital expandor with DLA (a = 0.5), mid- 
riser (c = 0.5), m = 3, and n = 4. 

and e 2 , and inhibits SB, shifting and the SR output. No adder is needed 

since Q = 0. 

The compression algorithm (corresponding to Algorithm 3) for the 
A-law with DLA is obtained by replacing expressions (48) and (49) by 

W(L, Y) e [vN, (1 + v)N) 



and 



V = 3[W(l, Y)] - vN. 



This algorithm leads to an A-law compressor circuit shown in Fig. 18. 
The overall performance of this compressor is similar to the n-\aw 
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circuit in Fig. 16. We will indicate only the differences between them. 
The number of cells in SR, is 7. There is no adder in the input since 
Q = 0. The gating pulse T L - V is obtained by the same circuit as in 
the expandor (Fig. 15). To facilitate the 0th segment output L = 0, 
a reset pulse t- 7 is added to BC. 

V. SUMMARY 

A unified formulation of segment type companding laws has been 
derived from which the detailed quantization structure of codecs can 
be explicitly represented. The decoder output levels and the coder 
decision levels are expressed in terms of standard forms and parameters 
a, b, and c. The standard form represents a chord-type piecewise linear 
approximation to smooth /z-law or A-law. The parameter a represents 
the segment edge effects, b relates to coding procedure, and the param- 
eter c represents the centering. These parameters are typically: DLA 
(a = 0.5 for b = ±0.5), RLA (a = for b = +0.5, or a = 1 for b = 
—0.5), mid-tread (c = 0) and mid-riser (c = 0.5). The effect of c is 
important in small signal range levels. The effect of a on the signal 
quality is very small, but the mixed use of DLA and RLA approaches 
introduce small degradations in S/N and net loss tracking. 

Synthesis of coders and decoders was described; and it was shown 
that the design figures such as amount of biasing in the decoder, the 



SR, 



CLOCK i 
INPUT Y(Z)' 



tt, a 2 <>3 «4 «5 a 6 °7 



t-7 



BC 



RESET DISABLE 
FOR t ~t 6 



RESET TO III 

DOWN 



T L -i 



V 



READ L AT tr 



r 





SR 2 




— ► 


a a 


"9 


a io 


'Ml 


T 

x 6 - 


1 ! ! 

e 5 e 6 e 7 y 



READ V AT t 6 



Fig. 18 — Thirteen-segment A-law digital compressor with DLA (o = 0.5), mid- 
riser (c = 0.5), to = 3, and n = 4. 
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type of logic to operate it, and the coding algorithm for the sequential 
comparison coder can be derived from our formulation. Bipolar codecs 
are also given in terms of the signed binary code representations. 

Another important consequence of this formulation is the systematic 
synthesis of digital compandors. The expansion algorithm is relatively 
straightforward. The compression algorithm for DLA was shown to 
be equivalent to employing the rule of "rounding" for the fractional 
digits of variable V. Examples of compandor circuits synthesized from 
these algorithms are shown for the segment /*-law and segment A-law. 

The A-law companding circuit appears to be slightly simpler than 
that for the ju-law, but the difference of the two is not so great as to 
weigh heavily in the choice between the /u-law and the A-law. It is 
noted that the formulation and the conversion algorithms developed 
here can also be conveniently applied for the synthesis of a digital 
converter between p-law and A-law or a digital attenuator. 
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APPENDIX A 

Proof of Congruence of Equations (15) and (16) 
For nonedge points X{L, V), from equations (7) and (14) obviously 

y d (X | b = +0.5) = y d (X + 1 \b = -0.5) (56) 

and, hence, (15) and (16) are congruent for nonedge points. For edge 
points 

y d (L, N-l\b= +0.5) = f(L, N + a- 0.5) -a + c, (57) 

y d (L + 1, | b = -0.5) = f(L + 1, a - 0.5) -a + c 

= f[L, N + 2(o - 0.5)] -a + c. (58) 

Here, these two equations are linear functions of a. Therefore, the two 
equations are equal if and only if a = 0.5. Thus, equations (15) and 
(16) are congruent for any X if and only if a = 0.5. 
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APPENDIX B 



Proof of the Statements in Section 2.S 

For given output levels the optimum decision levels lie halfway 
between two adjacent output levels. 10 At the nonedge points the proof 
is obvious from equation (14). At the edge points, the average of 
y (L, N - 1) and y (L + 1, 0) is given by 

\{y {L, N - 1) + y (L +1,0)} = y (L, N - 1) + * K(L) 

= y (L,N - J) + iA„(L)-(l + a) 



= yAL,N- 1+^ 



b = 4-0.5 



(59) 

This is equal to y d (L, N — 1 \b = 4-0.5) if and only if a = 0. Also 
equation (59) is expressed as 

y (L + 1, 0) - | X(L) = y (L + 1, 0) - \ A„(L + 1)-(1 + a) 

= y d [L+ 1,(1 - o)/4 | b = -0.5}. 

Therefore, this is equal to y d (L + 1, | b = —0.5) if and only if a = 1. 
This proves that the RLA condition is satisfied if and only if a = 
for 6 = +0.5 or a = 1 for b = -0.5. 

Next, for given decision levels the optimum output levels should 
be placed halfway between two adjacent decision levels. 10 At the 
nonedge points this is obviously satisfied. For edge points, since the 
half way point is equal to y(L. , V ,), from equation (23) it is shown 
that y = j/o if and only if a = 0.5. Therefore, the DLA condition is 
satisfied if and only if a = 0.5. 

appendix c 

Derivation of Algorithm 3 

Since in DLA (a = 0.5) the upper and lower bound approaches are 
identical, we will prove Algorithm 3 only for the case of b = —0.5. 
From Algorithm 2, equation (46) can be rewritten as 

y d (L, v) ^ y < y d (l, f + i). 

Since a = 0.5, this is valid for all V £ {0, 1, • • • , N — 1}. From equa- 
tions (56) and (27) the above inequality can be written as 

f + N g W(L, Y) < f + N + 1. (60) 
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Considering the range of t G {0, 1, ■ • • , N - 1}, we have as the L 
determining criterion 

AT S W(L, Y) < 2N. 
Then, from equations (50) and (60) we immediately have 

V = Z[Wil, Y)] - N. 
This completes the derivation of Algorithm 3. 

REFERENCES 

1. Smith, B., "Instantaneous Companding of Quantized Signals," B.S.T.J., 86, 

No. 3 (May 1957), pp. 653-709. 

2. Mann, H., Straube, H. M., and Villars, C. P., "A Companded Coder for an 

Experimental PCM Terminal," B.S.T.J., tf, No. 1 (January 1962), pp. 
1173—226 

3. Schreiner, S. M., and Vallarino, A. R., "48-Channel PCM System," 1957 IRE 

Nat. Conv. Rec, Part 8, pp. 141-149. 

4. Chatelon, A., "Application of Pulse Code Modulation to an Integrated Tele- 

phone Network, Part 2— Transmission and Encoding," ITT Elec. Commun., 
88, No. 1 (January 1963), pp. 32-43. .„.,,* 

5. CCITT (International Telegraph and Telephone Consultative Committee) 

Document COM XV-No. 77/E, "Interim Report of Working Party on 
Question 33/XV," 1966. 

6. Bucci, W., "PCM: A Global Scramble for Systems Compatibility," Electronics, 

42, No. 13 (June 23, 1969), pp. 94-102. 
7 Henning, H. H., "96 Channel PCM Channel Bank," 1969 Int. Commun. 
Conf., Denver, Colorado, June 9-11, 1969, pp. 34.17-34.22. 

8. Montgomery, W. L., "Digitally Linearizable Compandors with Comments on 

Project for a Digital Telephone Network," IEEE Trans, on Commun. 
Technology," Com-18, No. 1 (February 1970), pp. 1-4. 

9. Purton, R. F., "A Survey of Telephone Speech-Signal Statistics and Their 

Significance in the Choice of PCM Companding Law," Proc. IEE, 109, No. 1 
(January 1962), pp. 60-66. tt ^ 

10 Max, Joel. "Quantizing for Minimum Distortion," IRE Trans, on Inform. 
Theory, IT-6, No. 1 (March 1960), pp. 6-12. t 

11. Hoeschele, D. F. Jr., "Analog to Digital/Digital to Analog Conversion Tech- 

niques," New York: John Wiley, 1968. 

12. Chu, Y., "Digital Computer Design Fundamentals," New York: McGraw- 

13. Jackson, L. B.. Kaiser, J. F., and McDonald, H. S., "An Approach to the 

Implementation of Digital Filters," IEEE Trans, on Audio and Electro- 
acoustics, AU-16, No. 3 (September 1968), pp. 413-421. 

14. Kaneko, H., unpublished work. 

15. Schaefer, D. H. f "Logarithmic Compression of Binary Numbers, Froc. IKfc, 

49, No. 7 (July 1961), p. 1219. 



